Bethe-Salpeter equations for the collective-mode spectrum of a superfluid Fermi gas 

in a moving optical lattice 
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We derive the Bethe-Salpeter (BS) equations for the collective-mode spectrum of superfluid Fermi 
gases of equal mixture of atomic Fermi gas of two hyperfine states loaded into a moving optical 
lattice. In a moving lattice the superfluid state is unstable due to spontaneous emission of the short- 
wavelength rotonlike excitations which appear in the spectrum of the collective modes. It is shown 
that the spectrum obtained by the BS equations is in an excellent agreement with the collective- 
mode dispersion calculated by the perturbation approach, while there are some differences between 
the results obtained by density response function method and by the BS approach. The difference 
increases with increasing the lattice velocity, which can be seen in the analytical approximations for 
the dispersion relation in the long-wavelength limit in a weak-coupling regime. 
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I. INTRODUCTION 

The experimental ability to create a superfluid alkali 
atom Fermi gas in optical lattices^ opens a new op- 
portunity to study strongly correlated quantum many- 
particle systems and to emulate high-temperature super- 
conductors. Theoretical description of ultracold atomic 
Bose and Fermi gases in optical lattices has attracted 
much attention^r— Near the Feshbach resonance the 
atom-atom interaction can be manipulated in a control- 
lable way because the scattering length a s can be changed 
from the BCS side (negative values) to the BEC side 
(positive values) reaching very large values close to res- 
onance. On BEC side of the resonance the spin-up and 
spin-down atoms can form diatomic molecules, and these 
bosonic molecules can undergo a BEC at low enough 
temperature^— In what follow we focus our attention 
on the BCS side (negative scattering length) where BCS 
superfluidity is expected analogous to superconductivity. 
In particular, we consider an equal mixture of atomic 
Fermi gas of two hyperfine states with contact interaction 
loaded into an optical lattice. The two hyperfine states 
are described by pseudospins a =t,i- We also assume 
that the number of atoms in each hyperfine state per 
site (the filling factor) is smaller than unity, and that the 
lattice potential is sufficiently deep such that the tight- 
binding approximation is valid. The system in this case 
is well described by the single-band attractive Hubbard 
model: 



H 



(1) 



<l,J>,<T 



Here, the Fermi operator ipj a (ipi y(J ) creates (destroys) a 
fermion on the lattice site i with pseudospin er = j" , J, and 
Tt%,(j = tpj a ^pi,a is the density operator on site i with a 
position vector r^. fi is the chemical potential, and the 
symbol Yl<ij> m eans sum over nearest-neighbor sites. J 
is the tunneling strength of the atoms between nearest- 
neighbor sites, and U is the on-site interaction. On the 
BCS side the interaction parameter U is positive (the 



atomic interaction is attractive). We assume h — 1 and 
the lattice constant a = 1. 

In the case when the periodic array of microtraps is 
generated by counter propagating laser beams with dif- 
ferent frequencies, the lattice is not stationary anymore. 
It is expected that the formation of a BCS superfluidity 
is possible, but due to the presence of supercurrent (or 
quasimomentum q = tov, where to is the mass of the 
trapped atoms and v is the lattice velocity) the super- 
flow can break down. Recently, the stability of superfluid 
Fermi gases loaded into an optical lattice in the presence 
of supercurrent has been studied using the second-order 
time-dependent perturbation theory 2 ^ (this method is 
a generalization of the earlier approach by Belkhir and 
Randeria3^) , and the Green's function formalismi 24 i 25 In 
what follows, we derive the Bethe-Salpeter (BS) equa- 
tions for the spectrum of the collective modes in the case 
of a moving optical lattice, and compare the BS formal- 
ism with the above-mentioned two methods. 

This paper is organized as follows. In Sec. II, we derive 
the BS equations for the dispersion of the collective ex- 
citations. In Sec. Ill, the BS approach is compared with 
the results obtained by the perturbation theory, and by 
the Green's function formalism. 



II. BETHE-SALPETER EQUATIONS IN THE 
PRESENCE OF SUPERCURRENT 



We consider an imbalance mixture of 6 Li atomic Fermi 
gas of two hyperfine states \F = 1/2, m/ = ±1/2 >= 
\F = 1/2, my =fl> with contact interaction loaded into 
a moving optical lattice. The total number of 6 Li atoms 
is M, and they are distributed along N sites. The tight- 
binding form of the mean- field electron energy is £(k) = 
2 J (1 — cos k v d) — /i. In the presence of supercurrent 
with velocity v, the Cooper's pair has quasimomentum 
2tov. 

The Fourier transform of the single-particle Green's 



2 



function is a 2 x 2 matrix G = 



G t,t G t4 



In the 



G 4-4 Qi,i 

mean-field approximation the corresponding matrix ele- 
ments are as follows: 



G™(k,icj m ) = 
G^(k,iu m ) = 



u»(k) 




iu m -w+(k, q) 


zcj m + w_(k,q) 


<(k) 


u»(k) 


iu m -w+(k, q) 


iw TO + w_(k,q) 


G^(k, lWm ) 




1 


1 


_iw m - w+(k,q) 


iu m + w_(k, q) 



u q (k)v q (k) 



The symbol w m denotes w m = (2-k / (3){m + 1/2); m = 
0,±1,±2,..., /3 = {ksT)- 1 , k B is the Boltzmann 
constant, T is the temperature. As can be seen, 
the one-particle excitations in mean-field approximation 
are coherent combinations of electronlike w+ (k, q) = 
E(k,q) + r?(k, q) and holelike w_(k, q) = i?(k,q) — 
r;(k,q) excitations. The coherent factors it q (k) 
and v q (k) give the probability amplitudes of these 
states in the actual mixture. Here, E(k, q) = 



X(k,q) 
£(k,q) 



v/ X 2 (k,q) + A2, Uq (k) - 

and we have used the following no 



X(k,q) 
£(k,q) 



tations r?(k,q) = \ [£(k + q) - £(q - k)], x(k,q) = 
|E(q + k)+£(q-k)]. 

The number and the gap equations in this case are as 
follows (/ = M/N is the filling factor): 



1-/ 



N £-~t 



X(k,q) 



N^E(k,q) 



E 



i 



u 



7V^2£(k,q)' 



(2) 



The BS equations for the collective mode uj = w q (Q) 
and corresponding BS amplitudes 



* q (k,Q) 



/* q .t(k,Q)\ 
* q 4 (k,Q) 
*™(k,Q) 

\^(k,Q)J 



in the general random phase approximation (GRPA) are 
as follows: 



* q (k,Q) 



J2 *q(P, Q) + J M Y, $ i(P' Q)' 
P 



The first and second terms in the last equation represent 
the direct and exchange interactions: 
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( K^'^\k,Q,^ P ), ^ ,t4,t) (k,Q,tw I ,) 0\ 



K, 



4444) 
4444) 
4444) 



(k,Q,«jp), ^ q t,t44) (k,Q,ia;p) 



iC'" ; (k,Q, lWp ), X q 



4444) 



(k,Q,«j p ) 
Qc,Q,iu p ), X-^ ,t ' u) (k,Q,«j p ) 0/ 



,M 



/0 tf q ^Q,^), if q ^' t,t ' t) (k,Q,zu; p )\ 



if, 
K, 



4444)/ 



(k,Q, «*;„), if4 t ' t,t,J ' ) (k 1 Q,KJp) 



4444) 
q 

44-444) 



(k,Q,iw„), 4 t ' t,t ' t) (k,Q,za; p ) 
4444), 



\0 X-^" w " w (k,Q,«jp), if q + ' M '^(k,Q, lWp ) / 



Here, w p = (2ir//3)p;p — 0, ±1, ±2, ... is a Bose frequency, and we have introduced the two-particle propagator X^'''' k 
(i,j,k,l = {f,i})- 

A„"^;k. Q, iwt) = Y, G q J (k + Q; iu p + iu m )G% l (k; iuj m ) 



As in the case of a stationary lattice ) 20 i 21 the BS equations can be reduced to a set of two equations. At zero 
temperature, the corresponding equations are as follows: 



[w + n q (k,Q)-e q (k,Q)]G+(k,Q) 



U 
2N 



G+(P,Q) 



U 



G q (P, Q) - 7^7 E 7k,Q7^, Q (G+( P , Q) - G q (p, Q)) 



2N 



(3) 



4 E TO k,Q m p,Q [G q (P, Q) + G q (p, Q)] 



2N 



3 



[ w + tt q (k,Q) +£q (k,Q)] G"(k,Q) 
U 



u 

~2N 



Gq(P.Q) 



^ E [<q^,q - M, Q ] g+(p, Q) - ^ E tSqt5,q (p. Q) - G q ^ Q)) 



(4) 



^7 E <Q TO P,Q K (P. Q) + G q (P. Q)] 



2iV 



Here 



e q (k,Q) = £(k + Q,q) + £(k,q), 



f2 q (k, Q) = 2 [sin(fc„d) — sin(Q u d + k v d)] sm(q u d), 

V 

and the form factors are defined as follows: 

q q q i q q /q qq qq 

7k, Q — U k w k+Q W k u k+Q' J k,Q — M k M k+Q — ^k^k+Q' 



~q q q q x 

7k, Q — W k w k+Q — U k+Q u k' m k,Q 



q q q 



q q 

U k+Q v k- 



The BS equations ©- Q lead to a set of four coupled 
linear homogeneous equations. The existence of a non- 
trivial solution requires that the secular determinant 



u- 



7 q 

J q ~ 

7,7 

7 q 

7,771 



7 q 

7,7 



7 q ~ 

7 q 



7 q 

7,™ 



7 q ~ 

7,7 

7 q 7 q 

' 7,7 7 , 771 

J3 rj^ + xs 



(5) 

is equal to zero. Here we have introduced symbols 7 q b 
and J q ,: 



a/A 




FIG. 1: Collective mode energy in a stationary ID optical 
lattice at coupling U = 2 J. The filling factor, the gap and the 
chemical potential are / = 0.9, A = 0.346J and /i = 1.694J, 
respectively. The solid line is obtained by the Bethe-Salpeter 
method. The dotted line represents the spectrum obtained 
by the Green's function method (see FIG. 6a in Ref. [2^]) 
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'-4(k,Q)' 
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(6) 



•n q (k,Q)]' 



4(k,Q) 



The quantities q and 6^ q could be one of the four 
form factors: ; q q , m q Q , 7 q Q or 7 q Q . 

The secular determinant ([5]) is the main result in our 
BS approach. This determinant provides the spectrum of 
the collective excitations (in the particle and spin chan- 
nels) w(Q) in a uniform manner. 



Randeria— calculated the collective modes by lineariz- 
ing the Anderson- Rickayzen equations. In the GRPA 
the Anderson- Rickayzen equations are reduced to a 
set of three coupled equations and the collective-mode 
spectrum is obtained by solving the secular equation 
Det\A\ = 0, where 



A 



U 



Ij,l 



1,1 Jy,l 

u- 1 + I, 



J, 



7,m 



7,7 



J, 



u- 



7,m 



(7) 



III. COMPARISON WITH OTHER 
APPROACHES 

There exist another two methods that can be used 
to calculate the spectrum of the collective excita- 
tions of Hamiltonian (UJ). Decades ago, Belkhir and 



The quantities I and J are defined by Eq. when 
q = 0. The perturbation approach by Ganesh et al. 
provides the following 3x3 secular determinant (see Eqs. 
(B8) in Ref. [Hj], where there is a negative sign in front 



of the sum in the definition of Xq' 3 ) 



4 



^" 1 + (^+^ 7 q 7 -2J 7 q /)/ 2 (^t-^mV 2 • (8) 

- ^ (' 7 q 7 - 4)/ 2 ^ + ( J M + 7 7 q 7 + 2 ^V 2 / 



Since Det\D\ = det\A\, where 

u- 1 



Det\A\ = 



7 q 

w 7,m 



7,7 



7 q 



71 

7 ,m 

7 q 



IJ- 



m.in 



(9) 



one can say that the perturbation method is a general- 
ization of the Belkhir and Randeria method to the case 
of moving optical lattice. 



The secular determinant ([5]) can be rearranged as 



A B T 
B C 



0, where B T means the transpose matrix, 



and B and C are 1x3 and lxl blocks, respectively: 



B 



7 J^f,m I 1 



c 



(10) 



Thus, the BS secular equation assumes the form Det\A — 
B T C^ 1 B\ ~ 0. It should be mentioned that numerical 
calculations show that the contributions to A due to the 
term B T C~ X B are negligible, and therefore, there is an 
excellent agreement between the dispersions obtained by 
Ganesh et al. in Ref. [22j and by the BS secular equation 
not only in the case of a stationary lattice (see the con- 
clusion section in Ref. |H|), but in the case of a moving 
lattice as well. 

The method used by Yunomae et al.,— known as the 
Green's function approach, is a generalization of the Cote 
and Griffin^ approach to the collective excitations in 
s-wave layered superconductors. According to this ap- 
proach the spectrum of the collective excitations of an 
interacting electron gas (the model includes a Coulomb 
interaction v(r) — e 2 /|r| and a short-range attractive in- 
teraction g(r) between electrons) is obtained from the 
poles of the density and spin response functions in the 
GRPA. The response functions are obtained by using the 
Kadanoff and Baym formalism. The Kadanoff and Baym 
formalism is equivalent to the BS formalism. The BS 
equation for the two-particle Green's function is K~ 1 — 
j((0)-i _ wnere j((o) j s t ne two-particle propaga- 
tor, and the kernel I = I c i + I e xc is a sum of two parts. 
Id takes into account the direct interaction between the 
electrons, while I exc describes their exchange interaction. 
Cote and Griffin have ignored the long-range Coulomb 
interaction keeping in Id only the ladder diagrams in- 
volving the short-range interaction. The exchange in- 
teraction involves bubble diagrams with respect to both 
unscreened Coulomb interaction v(r) and short-range in- 
teraction attractive g(r). Thus, Eq. (2.33) in Ref. [26[ 
is equivalent to the BS equation K = + K^I d K, 





FIG. 2: Collective mode energy in a 2D optical lattice at 
coupling U — 4.5J and filling factor / = 0.8. The superfluid 
flows along (n, n) direction, and Q x — Q y — Q. (a) Stationary 
lattice (the gap and the chemical potential are A = 1.60J and 
(i — 3. 32 J), (b) Moving lattice q — 0.260 (the gap and the 
chemical potential are A = 1.60J and [i = 3.32J). Solid lines 
are obtained by the Bethe-Salpeter method while the dotted 
lines represent the spectrum obtained by the Green's function 
method (see FIG. 11 in Ref. [||). 



and Eq. (2.32) in Ref. [26| corresponds to the BS equa- 
tion K = K + KI exc K. In other words, the BS equation 
K _1 = K^^ 1 — Id — Iexc is equivalent to both, Eqs. 
(2.32) and (2.33) in the paper by Cote and Griffin. The 
kernel of the BS equations in the case when v(r) = and 
g(r) replaced by a contact interaction^ does take into 
account the direct interaction (terms which depend on 7 
and /) and exchange interaction (terms which depend on 
7 and m), and therefore, the Cote and Griffin results in 



5 



the case of an attractive Hubbard model should be sim- 
ilar to the results obtained by the BS approach. This 
statement is supported by FIG. 1 and FIG. 2. The two 
approaches provide very similar results, but the difference 
between them increases with increasing the lattice veloc- 
ity. To illustrate this statement, we compare the two an- 
alytical results for the speed of sound in a weak-coupling 
regime in a moving ID optical lattice. By following the 
calculations that have already been given by Belkhir and 
Randeria^ we reduce our secular determinant ([5]) to the 
following 2x2 determinant: 



2 q(^ 2 -4./^+A 2 

ma 

2A 



uj a ua 

4A 2 " 2A 

I- a 



where a — U/^y/ijjl — /i 2 . The solution w(q) = cQ 
provides the speed of sound which in a stationary lat- 
tice is c = vl — a^JVp — A 2 , where the Fermi velocity 
is Vf = (a/h)y/4Jfi — /i 2 . The same result has been ob- 
tained in Ref. 25]. In the case of a moving lattice, the 
corresponding determinant is: 



An 

uia , (2J—n)Qa tan(q) 
" 2A 2A 



" 2A 



(2 J— n)Qa tan(g) 
2A 

1 - a 



= 0. 
(12) 



0, (11) where a = U/(k^/U 2 cos 2 (<?) - (2 J - ^) 2 ) and 



An 



Q 2 a([i 2 - AJfi + A 2 ) uj 2 a (2 J - p,)Qauj tan(g) 1 2 2 [2J 2 - A + 4J 2 cos(2g)]Q 2 atan 2 (g) 



4A 2 cos 2 (g) 4A 1 



2A 5 



2A 2 



The last determinant provides the following expression for the phononlike dispersion in the long-wavelength limit 



(13) 



lj(Q) = (2J-/Z) tsm(q)Qa + QhVT~^ aJv$ + 4 J 2 [2 - 3 cos 2 (<?)] tan 2 (g) - A 2 



1 ~ a Or s 

1 H tan 2 (g) 



(14) 



Here, Vp = (a/h)y/ '4J 2 cos 2 (g) — (2 J — n) 2 is the Fermi velocity in the presence of a supercurrent. The speed of 
sound obtained by the Green's function method is)2£ 



oj(Q) = (2 J - pL) t&n{q)Qa + QhVl~a J - A 2 



1 - a 



tan 2 (g) 



(15) 



It should be mentioned that in the case when A/ J < 
1, the term 4 J 2 [2 — 3 cos 2 (g)] tan 2 (q) is more important 
than the term A 2 (l — a) tan 2 (q)/a in (|14l) and the term 
A 2 (l - a) tan 2 (q)/a 2 in (IT31) . 

IV. CONCLUSION 

We have derived the BS equations for the spectrum of 
the collective modes in the case when an equal mixture 
of atomic Fermi gas of two hyperfine states with contact 
interaction is loaded into a moving optical lattice. It is 



worth mentioning that there are two advantages of the 
BS formalism over the Green's function approach. First, 
within the BS approach we obtain the poles of density 
and spin response functions in a uniform manner, i.e. 
one secular equation provides not only the poles of the 
density response function, but the poles of the spin re- 
sponse function as well. Second, to obtain the poles of 
the density (spin) response function within the Green's 
function method, one has first to take into account the 
ladder diagrams with respect to the short-range interac- 
tion, and after that to sum up the bubble diagrams with 
respect to both, the short- and long-range interactions. 
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